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ABSTRACT 

Lyman-a absorption systems at z ~ 3 with iV HI ^3 x 10 14 cm -2 have been found 
to be enriched with a mean metallicity of Z/Zq ~ 10 -2 5 , and a large scatter in 
the metallicity. It is reasonable to assume that the process of initial enrichment of 
the intergalactic medium (IGM) at z ^ 3 also produced dust grains. We explore the 
implications of the presence of dust grains in the IGM at high redshift, in particular, 
the contribution of photoelectric emission from grains by hard background photons to 
the net heating rate of the IGM. 

We show that (i) the charge on dust particles and the characteristics of photo- 
emitted electrons differ substantially from those in the interstellar medium (ISM) in 
several respects: a) grains are exposed to and charged by photons beyond the Lyman 
limit, and b) because of this, the photoelectrons have typical energy of tens of eV. We 
also show that ii) silicates are more efficient heating agents than graphites; iii) small 
grains contribute mostly to the net heating; iv) at densities typical of the IGM at z ~ 3 
and for Lya absorbers, dust heating can be comparable or exceed photoionization 
heating within an order of magnitude; v) this increases the temperature of overdense 
regions, compared to the case of no dust heating, by factor of ~ 2. We discuss the 
implications of this extra heating source in Lya absorbing systems. 

Subject headings: cosmology : early universe - galaxies : intergalactic medium - 
galaxies : quasar : absorption lines 

1. Introduction 

Recent observations with high resolution have provided evidences for the existence of heavy 
elements in Lya forest clouds at redshifts z ~ 3 (Cowie et al., 1995, Tytler et al., 1995). Songaila 
and Cowie (1996) claimed that more than 75% of all Lya absorbers with N H1 ^ 10 14 cm -2 at 
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z ~ 3 are enriched. Tytler et al. (1995) reported similar results : all Lya systems with N ul £ 10 15 
cm~ 2 and 60% with iV HI ^ 10 14 cm -2 . These values are expected to be lower limits, since the 
measurements of heavy elements are S/N limited, and all Lya systems may contain heavy elements 
(see also Tytler 1988). 

The abundance of metals, however, remains somewhat uncertain, especially because it 
depends on the value of the ionization parameter (r, the ratio of ionizing photon density to 
the particle density) assumed. Nevertheless, for a reasonable range of T = l(T L5 -l(r 2 - 5 , the 
inferred metallicity is almost independent of the H I column density, and is of order Z/Z & ~ 1CT 2 
(Songaila & Cowie 1996). Rauch et al. (1997) found from the comparison of the data with 
their numerical simulations that the metallicity has a mean ~ 1CP 2 ' 5 solar and that there is 
a large scatter. Songaila (1997) has recently reported a lower limit to the metallicity of order 
Qmetais/^b ~ 1-4 x 10~ 5 h (1 + 2q z) 1 / 2 . For a low q Q universe (q = 0.02), h = 0.65 and a solar 
metallicity of 0.019, this implies a minimum metallicity of order 10 _3 3 -10 -2 ' 8 solar. 

Recent studies of metallicity and dust content in nearby galaxies have demonstrated that 
dust-to-gas ratio k and metallicity ( = Z/Z & correlate (Schmidt & Boiler, 1993, Lisenfeld & 
Ferrara, 1997). Though the correlation is not always good, it definitely shows that the production 
of metals in galaxies is accompanied with dust formation. Moreover, weak correlation between 
k and £ in blue compact dwarf galaxies may indicate that galaxies during active phases of star 
formation eject dust along with metals in gaseous outflows (Lisenfeld & Ferrara, 1997). It is, 
therefore, reasonable to expect that if the IGM and Lya absorbers are enriched with metals, they 
are also contaminated with dust grains, produced inside the galaxies and transported into the 
IGM along with metals. The presence of dust in the IGM and its consequences was first discussed 
by Rowan- Robinson et al. (1979), who suggested that dust in the early universe at z ~ 200 could 
lead to distortions of CMBR (cosmic microwave background radiation) (see also Vainer 1990; 
Loeb &: Haiman 1997). More recently, Fall et al. (1989) and Pei et al. (1991) have argued that 
the reddening of background quasars indicates directly that dust grains exist in damped Lya 
systems. Their dust-to-gas ratio ranges from k = 0.06 to k = 0.35 for different extinction curves 
assumed (here k = W 21 tb/N hi , kc = 0.79 for the Galaxy, Fall & Pei, 1989). Arguments based on 
the extinction of light by dust have led Zou et al. (1997) to similar estimates of k in damped Lya 
clouds toward the lensed quasar 0957+561. They explain the difference in the continuum shape 
of the two spectra by differential reddening due to the presence of dust grains. Recently, however, 
Pettini et al. (1997) have shown that the dust-to-metals ratio in damped Lya systems is only 
~ 1/2 of its value in local interstellar clouds (thus, dust-to-gas ratio is k ~ kcC/2). They infer 
this from the fact that Cr and other refractory elements are depleted only by a factor of 2 in 18 
studied damped Lya systems, much less than in the interstellar medium, and assuming that this 
depletion is due to freezing on dust grains. 

In this paper we assume that dust particles are present not only in damped Lya systems, but 
also in Lya forest clouds, following the idea that dust is associated with metals. We concentrate 
here on detailed study of dust charging in the high redshift IGM illuminated by the background 
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UV and X-ray photons and heating due to injected supra-thermal photo-electrons — in several 
respects these processes differ from those in the local interstellar medium. 

We do not discuss here how metals and dust particles are supplied to the IGM; we refer the 
reader to the papers where the IGM enrichment has been associated with population III stars 
(Ciardi & Ferrara, 1997, Gnedin & Ostriker 1997), early-forming very massive objects (Carr, 
Bond & Arnett 1984), galactic winds from galaxies at high redshifts (Silk, Wyse & Shields 1987; 
Miralda-Escude & Rees 1997; Nath & Trentham 1997). The estimates show that an enriched 
IGM with a metallicity of 10~ 2 ' 5 Z Q at z ~ 3 can be understood in terms of either Population III 
objects or galactic winds. 

The paper is structured as follows. In §2 we describe the basic processes and the equations, 
in §3 we describe the basic characteristics of charged high-redshift dust particles exposed to the 
background UV radiation, and photo-emitted electrons, and in §4 we present the results of the 
evolution of temperature in the IGM with photo-induced dust heating as an important energy 
source. We then briefly discuss the implication of our results for Lya absorbers. Throughout the 
paper we use a Hubble constant Hq = /150 50 km s _1 Mpc -1 and £7 = 1. 



2. The charge of a dust grain in the IGM 

2.1. The UV background radiation 

We first discuss the ambient photon background that is important for the process of charging 
of the dust at high redshift. The existence of a UV background radiation at high redshift has been 
inferred from proximity effect observations (Bajtlik et al. 1988). These observations show that the 
UV background radiation has an intensity at A = 912 A of Jgi2,-2i ~ 0.3-1, in the units of 10~ 21 
erg s- 1 cm.- 2 sr" 1 Hz" 1 (2 < z < 3) (Bechtold 1994; Cristiani et al. 1995; Cooke et al. 1996). 

On the theoretical front, following the pioneering work of Miralda-Escude &: Ostriker (1990), 
several authors have calculated the expected value of Jgi2,-2i from the quasars, using existing 
quasar luminosity functions, and taking into account the absorption of photons by the intervening 



Table 1 : Models for the UV background radiation 



A(A) 


QSO 


QSO+YG 


A = 912(A) 


J-21 = 0.5 


J-21 = 1. 


228 < A(A) < 912 


J v oc zv-°- 5 


J u oc u~ 2 - 5 


A = 228(A) 


J_2i = 0.025 


J_ 21 = 0.00123 


A < 228(A) 


J v oc v~ l 


J v oc v~ l 


Sl{= J912/J22&) 


20. 


813. 
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Lya absorption systems (e.g., Madau 1992; Haardt and Madau 1996). This predicted intensity, 
however, falls short of the inferred intensity from the proximity effect (Cooke, Espey and Carswell 
1996). It has, therefore, been suggested that there are other sources of UV photons at high 
redshift, e.g. young star forming galaxies, decaying neutrinos, or even a warm, collisionally ionized 
IGM. The predicted spectrum of the UV background in these cases however is much different, and, 
therefore, observational constraints on the spectrum hold the clue to the origin of the high redshift 
UV background radiation. The observations of He II absorption at high redshift put constraints 
on the spectrum of the ionizing background (Jakobsen et al.1994; Davidsen et al.1996). Meiksin 
& Madau (1994) showed that the He II opacity observed by Jakobsen et al.(1994) softness ratio, 
Sl = ^912/^228 ~ 40 at z ~ 3.3 assuming a reasonable value of the opacity due to diffuse H I (see 
also, Songaila, Hu & Cowie 1995). For the Davidsen et al.(1996) observations, Sethi and Nath 
(1997) found that the allowed range of models is still large, after taking various uncertainties into 
account (16 ^ Sl ^ 650). There has been recent claims that the spectral index changes around 
z ~ 3 (Songaila 1998, Reimers et al.1997, Hogan et al.1997, Savaglio et al.1997), but we assume a 
constant spectral index for simplicity. We discuss the effect of such a change in the section §3.3. 

For our purpose here, we will use two models for the UV background radiation for our 
calculation, bearing in mind the above uncertainties. The first one (referred to as QSO) represents 
the background radiation expected from only quasars, and the second one (QSO+YG), represents 
the radiation expected from quasars and young galaxies. The models are described in Table 1. 
The model QSO is a broken power-law fit to the spectrum obtained by Haardt and Madau (1996), 
for an intermediate redshift of z ~ 2.5. Note that their calculation was actually for q Q = 0.1. The 
model QSO+YG assumes the QSO luminosity function of Boyle (1991), as in Madau (1992), and a 
population of young galaxies, which is used to make Jgi.2,-21 = 1 using the A3 model of absorption 
of Giroux & Shapiro (1996) (again at z ~ 2.5; see also Sethi and Nath 1997). The values of the 
softness ratio Sl for each model is given in Table 1. Although the broken power laws are crude 
approximations to the actual spectrum, they suffice for our purpose of illustrating the effect of 
dust charging and heating at high redshift, in the light of the uncertainty in our knowledge of the 
UV background radiation. 

Beyond 228 A, we use J v oc u^ 1 till hv = 3 keV for both models. This again mimics well the 
Haardt & Madau (1996) spectra for 10 A ;$ A ^ 228 A. For ^ 3 keV, we use the quasar X-ray 
spectral index of 0.9. Zdziarski et al.(1995) found a good fit to the data of Seyferts with a spectral 
index of 0.9 till ~ 100 keV. However, photons with energies much above 1 keV do not contribute 
significantly in the charging and heating processes due to small absorption coefficients at these 
energies. 

The evolution of the background radiation in redshift also remains uncertain. For simplicity, 
we assume that Jgi.2,-21 remains a constant for 2 < z < 5, and Jgi.2,-21 oc (1 + z) 2 - 8 for z < 2. This 
is a fit to the evolution of Jgi.2,-21 as calculated by Madau (1992) for z < 2. Also, we assume that 
the spectral shape remains a constant in time, although in reality the spectrum becomes softer at 
higher redshift due to increase in absorption (Miralda-Escude & & Ostriker 1990). Again, however, 
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the uncertainties involved are large, and we assume a constant spectral shape for simplicity. 



2.2. Basic processes for charging grains 

In this paper we follow Draine (1978) and Draine & Salpeter (1979, hereafter DS) in describing 
the processes determining the grain charge. The equation governing the charge on a grain is 
written as follows 



dq 
~dt 



= 4ira eHjJj 



(1) 



where Jj is the partial current due to particle of the i-th kind, and a is the grain radius. The 
currents which mainly determine grain charge in the ISM are due to: i) impinging electrons, 
ii) impinging ions (mainly protons), Hi) UV photons. In principle, secondary processes can be 
associated with these three (such as secondary electron emission), however, in the conditions we 
are interested in they play a minor role. In case of incident electrons and protons, currents Jj are 
enhanced due to Coulomb interaction by factor g{x). For electrons, with x = eU/kT, where U is 
the grain potential, (Spitzer, 1978), 



9e(x) 



exp(x) if x < 
1 + x if x > . 



(2) 



For ions, it is given by, 



where z is the ion charge. 



9i{x) 




if x > 
if x < , 



(3) 



2.2.1. The current due to impinging electrons 

Grain particles in low density IGM exposed by UV and X-ray photons are positively charged, 
and the electron current can be written as (DS), 

/ kT \ 1 / 2 

Je = ~ ne {2^n~J «*e>-W)(l + *), (4) 

here (s e ), the thermally averaged sticking probability for incident electrons. This probability, 
(s e ) ~ 1 for thermal electrons with temperature T and grain size a satisfying the condition 
(Draine, 1989), 
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1.5 



" >10 {wk) A - (5) 

which is valid for temperature T < 10 5 K, the range of temperatures relevant in our case. (5) is 
the secondary emission coefficient, for incident thermal electrons with kT <C 300 — 500 eV (i5) < 1. 



2.2.2. The current due to impinging ions (protons) 
DS give, for the proton current, 

/ kT \ 1/2 

where (s p ) — 1, the effective sticking coefficient, (u p ) <C 1, the secondary electron emission 
coefficient, and f3(Y z ) <C 1, the charging due to sputtered charge particles, as can be estimated 
from DS for temperature T < 10 5 K. It is seen clearly that in fully ionized IGM, contribution of 
protons to grain charging is of order \J p \ ~ 0.03| J e \. We neglect the current from helium ions, since 
for the situations we are concerned with in this paper, the contribution from ions is in general 
small compared to electron and photoinduced currents. 



2.2.3. Photoelectric current per area 
Photoelectric current per grain area is 

oo 

Juv = J Q{a,v)y v ^-j-dv , (7) 

where J v is the spectral energy flux of incident photons ([erg cm~ 2 s _1 Hz -1 sr -1 ]), y v , the 
photo-yield, i.e. the number of photo-emitted electrons per absorbed photon, Q„{a), the efficiency 
of photoabsorptions, and v m { n = (w + eU) /h corresponds to the minimum energy for which 
electrons can escape from the surface of a dust grain, w being the work function. The normalized 
photo-yield can be written as (Draine, 1978, Bakes & Tielens, 1994, hereafter BT) , 



where (Draine, 1978), 



Vv = Vo{l-^)fz(a.), (8) 



h{a) = [a) (C 2 -2 C + 2- 2e-C) ' (9) 
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with 



a a 




(10) 



= T + P 



where a is the grain radius, l a , the photon attenuation length, l e , the electron escape length; BT 



the ionization potential of a grain I z ) value of yo is determined by Draine (1978) as yo = 0.5 to 
reproduce the enhanced value of the photo-yield of a small particle over the bulk graphite, and by 
BT as yo = 0.14 to fit experimental data for coronene (Verstraete et al, 1990). I z is the ionization 
potential for a grain with charge Z (BT), and is given by, 



Here e is in esu units. It is important to note here that when the photo-induced current is 
calculated for the ISM in the Galaxy, the upper limit in the integral is defined as the Lyman 
limit frequency ul (Drain, 1978, BT), regarding the fact that Lyman continuum photons are 
confined within HII regions. This restricts from above the charge of dust grains, the energy of 
photo-emitted electrons and thus the photoelectric heating. In our case, charging of dust grains 
and photoelectric heating are determined by photons with energy above Lyman limit. BT stress 
that near threshold the photo-yield for graphite grains is usually overestimated. In our case, with 
photon energies and electric potentials much higher than in the ISM, the contribution from the 
threshold region is small. 



For silicate particles Q can be reasonably approximated for sizes what we are concerning 
a < 0.1 /xm, and photon energies e > 10 eV by (Ferrara & Dettmar, 1994) 



give l a = 10 6 cm, and l e = 10 



cm. The asymptotic (at photon energies much larger than 



(ii) 



2.2.4- Photoab sorption efficiency 




(12) 



where 



7a = 2.5 (10a) ai - 1 



(13) 



here and throughout the paper a is in /im, if not specified. 



For graphite particles, we use the approximation, 



Q(a, e) = 30 x a exp(-ll. 9\fa)a\ (e) 
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+ 60 x ° exp(-7.5^)^2(e) 

10a . . . . 

+ (l + ^d + lOa) ^' (14) 

where <7j, the cross sections per carbon atom in a particle (in MegaBarns), are described below. 
For <7i and 02, we used the data given in BT (see their Fig. 1), which we approximated as 

, . _ f 5.24 x 10- 4 e 3 - 78 , e<17.29eV . . 

ai[€) ~ 1 1.37 x 10 5 e- 3 - 02 , e > 17.29 eV 5 { j 

, , f0.032e 2 ' 5 , e<6.13eV , 1C , 
a2(e) = { 112.7c-, £ >6.13eV ; (16) 

^(e) = 1.5 + 2 x 10~ 4 (e - 0.22)(e - 18)(e - 40) , (17) 

here e is given in eV. This approximation for Q(a,e) agrees with the results of Draine and Lee 
(1984) (see their Fig. 4) for grain sizes a = 0.03 — 0.1 and with the results of BT (see their Fig. 1), 
which are based on an extrapolation of the results of Draine and Lee to smaller particles, to 
within 10-15 %, for 3eV < e < 40 eV, and 20-30 % for e < 3eV. BT showed that the ionization 
cross sections obtained in this way are in good agreement with the data for small PAH coronene 
(Verstraete et al. 1990). 

At energies e > 40 eV, for Q(a, e), we use the approximation given by DS 

/40 eV\ 

Q(a, e) = Q{a, 40 eV) ; (18) 



1 keV \ 3 

Q(a,€)=Q(o,lkeV)(— — ) . (19) 



and for e > 1 keV, 

, ,. kc-V _ 

These approximations give reasonable fits (within 10-15 %) in the entire energy range. One 
should stress however, that in the high-energy interval, e £ 1 keV, the Q(a,e) oc e~ 3 dependence 
is contaminated by very sharp increases in Q when the photon energy approaches the edge of an 
atomic subshell, and can exceed the Q(a, e) oc e -3 value by 1-1.5 orders of magnitude, (Dwek & 
Smith, 1996). From this point of view our results for dust charging and heating can be considered 
conservative. 



2.3. Grain potential 

It is readily seen from the charge equation that characteristic charging times caused by the 
photo- induced and electron currents, for n e > 10~ 7 cm~ 3 , are much shorter than the Hubble time 
at z ;$ 20 (iff ~ 10 15 s) and sound crossing time for a cloud with R > 1 kpc with temperature 
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T ~ 10 4 — 10 5 K (tp, ~ 10 15 s). Therefore, the grain charge is kept at equilibrium, q = 0, and can 
be easily determined from 

n e ( (1 + x = / Q(a,v)y v —^dv + nJ . 20 

It is clear that the proton current can be comparable to the photo-induced current when the 
particle flux for protons ~ n p v p is comparable with the photon flux ~ J v /h. More precisely, for a 
power law spectrum J u oc v~ a it is equivalent to the condition, 



' SkT\ l/2 4irQyJ (hv L 



irm p I ah \ w 



(21) 



where Jo is the spectral energy flux at the Lyman limit frequency, vl, Q and y, the absorption 
efficiency and photo-yield averaged over the spectrum. For typical values Jo ~ 10~ 21 , Q ~ 0.1, 

— 1/2 

y ~ 0.1 this is equivalent to n e ~ 0.2T 4 . Thus, for conditions in the IGM at z ~ 3 photo- induced 
current dominates dust charging. 



2.3.1. Analytical estimates 

For grains exposed to hard UV radiation, the charge and potential are expected to be high, 
x S> 1, and can be estimated for a power-law spectrum of incident photons J oc v~ a and a 
power-law dependence for Q(a, v) oc Q{a,vi){y /vl)~^ \ as 



x ~ 



10 a+/3-2 

a+/3+1/2 Q(o,i/ £ ) 
n e T 4 '(« + /?) 



1/(q+/3+1) 

(22) 



Here, we have used y = 0.1; our assumption of x S> 1 is valid for n e < 10~ 3 — 10 -4 cm~ 3 , and 
T4 < 10. For a ~ 1, (3 ~ 3 (such /3 is typical for graphites particles of small sizes in the energy 
interval e = 17.3 — 40 eV), this gives a weak dependence of the potential on gas density and 
absorption characteristics, x ~ 2n~°' 2 T 4 -0 ' 9 [J_2i(3(a, i^l)] ' 2 . In Lya clouds with n e ~ 10~ 4 cm~ 3 
and T ~ 10 4 K x ~ 10[J_2iQ(a, ^l)] ' 2 which is much larger than the value for grains in the local 
ISM illuminated by photons below Lyman limit. For graphite particles with larger radius, a ~ 0.1, 
P ~ 1 - 1.5, and x ~ 0.8n-°- 25 r 4 -°- 6 [J_ 21 Q(a, ^)] " 25 - n- - 28 T 4 - a86 [J_ 21 g(a, ^ L )]°- 28 . 

Silicate particles have less strong dependence of Q(a, e) on photon energy in the interval 
e > 10 eV, (3 ~ — 1 for a = 0.1 — 0.001, respectively, which results in a stronger dependence 
of grain potential on the parameters of ambient medium: x ~ 0.3n~ a5 T 4 " a75 [J_2i<2(a, i^l)] ' 5 - 
0.8n~ a33 T 4 ~ a83 [J_2i<2(a, ^l)] ' 33 - Note that for the dependence of the graphites grain potential 
on Q(a) and on grain size is weaker than for silicates. 
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2.3.2. Numerical results 



Figures la and lb show the dependence of grain potential x = ell/kT for graphites and 
silicates, respectively, on the ambient density for hard (QSO, solid lines) and soft (QSO+YG, 
dashed lines) spectra of background photons, for a grain size of a = 10 -5 cm and for fixed gas 
temperature T = 10 4 K. Figures 2a and 2b show the dependence of x for graphites and silicates 
on the grain size a for different ambient densities for the QSO and QSO+YG spectra. 

The figures show that the dependences follow qualitatively the above simple estimates. In 
general, hard spectrum (QSO) results in higher value of x, especially for silicates, and both 
graphites and silicates have comparatively lower grain potentials for soft spectrum. QSO spectra 
with low equivalent a result in more strong dependence of grain potential on density than 
QSO+YG spectra with larger equivalent a do. Graphites do not show much dependence of the 
potential on the grain size, while silicates do, which is due to more weak dependence of Q(a, e) on 
photon energy for silicates, as seen from analytical estimates. It is also seen that grain potential 
on silicate particles depends on the ambient density stronger than on graphite particles. 

2.4. Photo-electric heating 

The heating rate due to photo-electric effect in dust grains (in erg cm -3 s _1 ) is given by, 



H p h = nil J daAf g (a)a(a) J dvQ{a,v)y v - - (hv — w — eU (a)) , (23) 

"ml(«) 

where v m i(a) = (w + eU + 3kT/2)/h, the minimum frequency of photons needed to eject electrons 
which have sufficient energy to heat the surrounding gas, N g {a) = dn/da, the size distribution of 
grains, is taken a power-law M g (a) = Aa~ q . The normalizing factor A can be determined from 



/4 
a~ q -TT Pg a 3 da = 10~ 2 Cm#n , (24) 



where a m i n = 5 x 10 7 cm, a max = 4 x 10 5 cm (see, O'Donnel & Mathis, 1997), are the minimal 
and maximal sizes of dust particles, respectively. Then, the coefficient A is 



A = 3.8x 10~ 27 



(4-9) ( 



Z 



ii 



4-g 
Ctmax 



4-q 



(25) 



for typical grain density p g = 1 g cm 3 . For U(a) in Eq. 23, we use the results for grain charge as 
obtained in 52. 3. 
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2. 4-1- Mean energy of photoelectrons 

An important parameter which determines the photo-electric heating is the mean energy of 
photo-emitted electrons, 

oo 

/ di/Q(a,i/)y v ^(hu-w-eU(a)) 

E(a) = ^ - . (26) 

/ dvQ(a,v)y v ^ 

In the local interstellar medium E(a) is of the order of hv^ — w ~ 5.6 eV. In a medium exposed 
to hard UV photons it is restricted only from below, and the exact value depends on the UV 
spectrum. Figures 3 and 4 show this dependence for graphite and silicate particles. It is seen that 
E(a) can reach hundreds of eV at low densities, while typical values for moderate densities is 
30 - 50 kT ~ tens of eV. 

Note that although the curves in Figures 3 & 4 show that in general E(a) for grains with 
higher charges is large, it should not be taken to mean that higher charge causes the mean electron 
energies to be larger. The dependence of E(a) on ell/kT at T = 10 4 K is shown in Figure 5 for 
silicates (solid lines) and graphites (dotted lines), for QSO type (thick lines) and QSO+YG type 
spectra (thin lines) (the upper set of curves are for a = 10~ 5 cm and the lower set of curves are 
for a = 10 -6 cm). It is seen that E(a) is a monotonically increasing function of eU /kT, and the 
nature of the curve depends on Q(a, v), y u and the photon spectrum. 

The differential heating rate (per unit grain) is a product of absorbed photons and E(a) 

oo 

^ = E(a) J d,Q(a,u)y v a(a)^ , 

and can be therefore roughly estimated as dH p h(a)/dn g (a) oc E(a)(eU) l3+a ~ 1 oc 
E{a)nf +a - 1)/[ P +a+1) , i.e. natter than linear. 



2.4-2. Differential heating rates 

Using the approximate solution for dust potential what we found in section 2.3.1 one can 
estimate partial contributions to the total heating rate from silicates and graphites. In the integral 

oo 

dH vh f , ~, s , AnJ u 



dn g {a) 

"ml(o) 



/47T J 
dvQ(a, u)y u a(a) - - (hi/ — w — eU(a)) , (27) 
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we can neglect w in comparison with hv and ell, and using as above power laws for the flux J and 
absorption coefficient Q(a, v) we get 



? ^ _4x 10- 5 /(a,/5V(«K^ +1 r 4 2(a+3+1, [«/-2iQ(a,^)]™, (28) 
an g {a) 

where f(a,f3) ~ 1. It is seen from this equation that the differential heating rate, both for 
graphite and silicate particles, depends on grain radius through Q{a,v£) rather weakly, and thus 
contribution from grains of small radius dominates due to the factor N g {a) oc a~ q . However, 
the integrated heating rate is insensitive to the grain size distribution: for Q 2 / a+ P +l ~ a 7 ?, with 
< r] < 1 the dependence of heating rate on q is close to be eliminated with normalizing factor A 



in Eq. 25 



In figure 6 we show heating rates (per unit volume) due to graphite and silicate particles, 
respectively, integrated over size distributions for both incident photon spectra, and for different 
size distribution of grains, q = 3.5 and q = 3.9. For comparison we show also the photoionization 
heating rate for a primordial gas. It is seen that photoelectric dust heating depends on density less 
strongly than photoionization heating. Silicates produce weaker dependence of heating rate on 
density, which can be explained by the weak dependence of grain potential and the mean energy of 
photo-electrons on density (see, Figs 2 and 3). Due to the fact that both ell and E(a) for silicates 
are higher than for graphites, the heating rate by silicate grains is also higher. Changing the size 
distribution results in a weak change (within factor of 1.5) of the heating rate: dust with q = 3.9 
heats more efficiently than with q = 3.5, as expected. In Figure 7, we show the heating rate 
as a function of the lower limit of the grain size assumed, for different grain types and ambient 
densities. The curves show that the heating is dominated by small grains in general. In brief, 
in all cases the photo-electric dust heating is comparable or even exceeds (at low densities) the 
photoionization heating. 

The dependence of the dust heating rate on the ionizing radiation flux at Lyman edge 
J{y£) is quite weak, particularly for QSO+YG spectra with large a. However, it can lead to 
variations in temperature due to local inhomogeneity in exposing UV flux. On the other hand, the 
photoionization heating does not depend on the flux in a fully ionized medium with njji <C njj 
for low density (n <C 10~ 35 cm -3 as considered here), as the photoionization heating rate is 
proportional to the recombination rate and thus weakly depends on the UV flux.. 

Dust heating can also stimulate at appropriate conditions spatial variations in temperature 
via thermal instability. In the simplest case of steady energy balance: A(T)n 2 = H p h, with H p h 



depending on density and temperature as given in Eq. |28|, Field's criterion for thermal instability 
of condensation mode (Field, 1965) is written as 

rflogA 3 a + /?-l 

dl^T- 1+ 2a + (3 + l <0 > (29) 
which shows that for a + (3 < 1 (hard spectrum and dominance of silicate grains with a > 0.01 /urn), 
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dust heating works as a destabilizing factor. 

3. Thermal evolution of the clumpy IGM with dust heating 

Next, we compute the thermal evolution of the clumps in the IGM with overdensities of order 
~ 10, which either keep their overdensity constant (i.e. , adiabatically expand), or are confined by 
some mechanism and keep their density constant. These are described below in detail as Models I 
and II. 

We assume that the gas in Lya systems evolves passively, driven by the UV background 
flux and the expansion (or non-expansion) of the cloud. This is consistent with the picture of 
Lya systems as described by the recent numerical simulations (e.g., Miralda-Escude et al.1996) 
in which Lya system gas is not involved in star formation and are not affected by the processes 
related to it. 

We note here that Miralda-Escude & Rees (1994) also calculated the evolution of temperature 
in overdense regions of IGM, which first collapse under gravity, during which the temperature 
rises to ~ 10 5 K due to adiabatic heating, and which then evolve with constant (proper) density, 
during which, due to high density and the consequent high recombination and line cooling rate, 
temperature falls down quickly to ~ 10 4 K.They, however, considered much higher densities 
( ~ 1CT 3 cm~ 3 ) than we use here ( ^ 1CT 4 cm~ 3 ) (at lower densities, the balance between the 
heating and cooling processes will be different, for example, adiabatic cooling will be more 
important). 

The evolution of temperature in an overdense region in the IGM would depend on the initial 
temperature (and epoch) assumed. If the gas in the overdense region were initially compressed 
adiabatically, then the maximum temperature would be T c ; ~ (pel/ Pigm) 2 ^ Tigm, where p c i 
is the density in the overdense region. For overdensities of order 5-10, that are considered to 
be typical values for Lya absorbing systems (see below), T c i ~ 3-5 Tigm- The value of Tigm, 
however, depends on the thermal history of the IGM (Miralda-Escude & Rees 1994) which remains 
uncertain. We, therefore, for the sake of illustrating the importance of dust heating, assume a few 
values of the initial temperatures of the overdense region Tj„ which are of order 10 4 K, and follow 
the thermal evolution in time. 



3.1. Models: heating and cooling processes 

In numerical calculations we used the heating rate due to photo-electric effect in dust grains 
as described above with grain size distribution M g oc a~ q with q = 3.5. For the relative fraction 
of graphites and silicates, we use the one given by Pei (1992) for the Large Magellanic Cloud 
(LMC), which seems suitable for low metallicity objects at high redshift (Pei, Fall & Bechtold 
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1991). We also show the thermal evolutions with only graphites and silicates involved. The case of 
only silicates correspond to the Small Magellanic Cloud (SMC) type dust grains (Pei 1992). For 
other heating processes apart from dust heating, we use photoionization heating of a primordial 
gas, for both models of UV radiation. The heating due to photoionization of carbon atoms is 
smaller than this by a factor of 3 (10) for the QSO+YG (QSO) spectrum for a metallicity of 
Z/Zq = 0.01. For cooling, we use recombination cooling, cooling due to collisional ionization, line 
cooling (in the presence of the UV radiation) and bremsstrahlung for a primordial gas, as given in 
Cen (1992). Cooling due to heavy elements at this metallicity and temperature is unimportant. 
Inverse Compton cooling at z ^ 5 is not important and is neglected here (the corresponding 
cooling time-scale is tic ~ 9 x 10 12 (1 + z)~ A yr; Tegmark, Silk & Evrard (1993)). 

In the Model I, we assume that the overdense regions adiabatically expand keeping the 
overdensity a constant in time. In this case, therefore, we also include cooling due to adiabatic 
expansion. In the Model II, we consider regions with constant (proper) density, in which the 
gas is confined by gravity or shock ram pressure, or some other confinement mechanism. For 
Model I, we use Pd/^bPc = 10, for all redshifts, where p c \ is the density inside the cloud, p c is the 
critical density and 0^ is the global mean baryon density in the units of p c . For Model II, we use 
p c i/^lbp c = 10 (at z = Zi n ) for different values of the initial redshift Zi n . We use = 0.03 and 
initial temperature Tj„ = 10 4 ,2 x 10 4 ,5 x 10 4 K. Our choice of the parameters follow the result 
of Cen & Simcoe (1997), that Lya absorbers with N H1 ~ 10 14 cm~ 2 had p c i/pb ~ 10 at z ~ 3 for 
Qb ~ 0.03. The minimum width of Lya lines at z ~ 4 has been recently found to be b = 15 km 
s _1 (Kim et al.1997). If the minimum width (in the distribution of b values) is considered to be 
due to gas temperature (see below for more discussion on this), this corresponds to a temperature 
of ~ 10 4 K and proves the presence of overdense regions with such temperatures at z ~ 4. Here 
we consider the temperature evolution of such regions in time. 

In reality, the gas in Lya systems is expected to behave in a way which lies between our 
two models. Numerical simulations show that the structures in the IGM corresponding to Lya 
forest absorption lines are sheet-like structures (Miralda-Escude et al.1996). The gas is confined 
in one dimension and flows along the sheet toward semi-spherical vortices. Therefore, the density 
is expected to fall slower than (1 + z) 3 . Also, since these structures are transient, lasting about a 
Hubble time, our calculations would be relevant only on time scales shorter than this. 

3.2. Models: numerical results 

We show the evolution of the gas temperature in Figure 8a and 8b for graphites and silicates. 
The dotted and solid lines refer to the cases without dust, and with dust with £ = Z/Zq = 0.01. 
The plots show the evolution for various initial temperatures, for the two models of background 
photons (upper panels for QSO and lower panels for QSO+YG spectra). The case of silicates, 
of course, correspond to the grain mixture of the SMC type (Pei 1992). In figure 9 we show the 
case of grain mixture of the LMC type (Pei 1992) for the QSO spectrum for three values of initial 
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temperatures and metallicities. 

It is seen that the dominant contribution to the heating comes from silicate grains for both 
type of spectra. This is due to the fact that for silicates, the absorption efficiency does not decrease 
with increasing photon energy as sharply as it does for graphites, and therefore, the high energy 
photons contribute substantially in this case. In case with hard background spectrum (QSO), the 
dust heating is more efficient for the same reasons — due to contribution of high energy photons, 
and as a result of the high mean energy of photo-emitted electrons. In Models I without dust 
heating adiabatic cooling decreases temperature sharply, so that at z = 1 — 2 it is below 10 4 K. 
Dust heating is able to support temperature at the level T ^ 10 4 down to z = 0. In Models II 
with dust heating temperature is factor of 1.5-2 higher than without dust heating. At z < 1 both 
photoionization and photoelectric heating decrease due to decrease of the UV flux. 

Since dust heating manifests itself better (compared to other sources of heating) at lower 
densities, it is conceivable that the diffuse, low density, 'intercloud' component of the IGM, if any, 
would have a much different thermal history with dust than with only photoionization heating. 
The density of such a component remains, however, uncertain. Recent numerical studies show 
that most of the baryons in the IGM are in the overdense regions, leaving about 5% in the 
diffuse component (Zhang et al.1997). Their results show a fragmented IGM at high redshift, 
and they claim that even the 'diffuse' medium is also possibly fragmented into scales smaller 
than the resolution of the simulation. For fif, = 0.03 and /150 = 1, this corresponds to a density 
~ 4 x 10~ 9 (1 + z) 3 cm~ 3 (using 5% of the baryons in the diffuse component). Judging from the 
heating curves of Fig. 6, such a component would be hotter than if only photoionization is used. 
However, the temperature would depend on the thermal history of the IGM. 

3.3. Uncertainties 

Here we discuss briefly the effect of uncertainties in the parameters on our results. Since 
the type of dust particles in Lya systems is uncertain, we showed the temperature evolution for 
different types of dust (graphite and silicate) in Figures 8a and 8b. It is seen that the effect 
on the temperature of the gas is by factors of order ^ 1.5. We also show the effect of different 
metallicities on the temperature evolution in Figure 9, for Z/Z & = 0.003,0.01. The difference 
in the temperature of the gas is again of order ~ 1.5. Since the actual temperature evolution 
of the gas would depend on the initial conditions, the readers are urged to use the heating rate 
shown in Figure 6, and scale it for the required metallicity for different cases (the heating rate is 
proportional to the metallicity). 

Finally, the evolution of the UV background radiation is still uncertain. For example, it 
has been claimed that the spectrum changes at z ~ 3 from soft at higher redshifts to a harder 
spectrum later (Songaila 1998 and the references therein). We have already shown the heating 
rates (Figure 6) and temperature evolution of the gas (Figures 8a and 8b) for a hard (QSO) and a 



-16- 



soft (QSO+YG) spectra, with the spectral index being constant in time. Figures 8a and 8b show- 
that the temperature of the gas at later redshifts (z ~ 2) differ by very small factors for different 
values of Zj n . Since the heating rate for the soft spectrum case is very small compared to the hard 
spectrum case, a change in the spectrum at z ~ 3 would essentially be equivalent to the case of 
Zm ~ 3, and as explained above, it would effect the temperature evolution at later redshifts only 
negligibly. For Zj n = 4, the temperature would evolve much like the case of just photoionization 
between 3 4 (dotted lines in Figures 8a and 8b). 

Similarly, the possibility that the value of J%n-i\ drops by a factor of 3-10 from z ~ 2.5 to 
z ~ 4.2 (Williger et al.1994, Kim et al.1997), would essentially mean a very small heating rate 
at z £ 3, but as explained above, it would not effect the temperature evolution at later redshifts 
much. 



4. Discussions 

Finally, we discuss various implications of our results on photoelectric heating due to dust in 
IGM. To begin with, we discuss the plausibility of having dust in the IGM in the first place, i.e. , 
the survival of dust grains in IGM. 

4.1. Survival of dust particles 

Before the discussion of the heating process due to dust, let us first consider the question of 
the survival of dust grains in high-redshift Lyman-a clouds. As already mentioned, these clouds 
have typical densities ^ 10~ 4 cm -3 at high redshift, metallicities ~ 1O~ 2 ' 5 Z , and they are 
exposed to background radiation ~ 10 -21 erg cm" 2 s _1 sr _1 . It is not clear if these clouds were 
enriched by metals driven out of earlier generation dwarf galaxies with galactic winds or these 
metals come from in situ star formation; though the small scatter in observed metallicities as iV HI 
changes by over four orders seems to favour the former explanation (Cowie & Songaila 1998). 

In the model of enriching the IGM by galactic winds at high redshifts (Silk, Wyse & Shields 
1987; Miralda-Escude & Rees 1997; Nath & Trentham 1997), the enriched gas is swept along 
with the ambient IGM by winds erupting out of dwarf galaxies due to the injection of energy by 
supernovae from an early bout of star formation. Dust grains could form either the stars or in the 
shells of supernovae (Dwek 1998). Inside the parent galaxies, the processes of destruction of dust 
grains could be as efficient as in the Galactic ISM. The destruction of dust grains in our Galaxy 
is still not well understood (see McKee et al. 1987, and references therein). However, detailed 
modelling show that it is reasonable to assume that a fraction (~ 10%) of refractory dust survives 
the destruction process in ISM (Dwek 1998). 

When the enriched material, presumably along with the dust grains, are ejected in the IGM, 



-17- 



the destruction processes become less important. For sputtering, for example, the grain life time 
(above the sputtering threshold) is, 

t sp ~ ^ , . 1 n 6xl0 13 ^s. (30) 

Y sp (v)n s a g Q coU (y sp /0.01) n s T 4 1/2 

Here, N a is the number of heavy atoms in the grain, Y sp < 0.01 is the sputtering yield (number of 
atoms ejected from the grain per impinging atom) (Draine & Salpeter 1979), < v > is the thermal 
speed of the ambient particles with density n s and temperature T = T4 10 4 K, a g = ira 2 is the 
geometrical cross-section of the grains with radius a and Q co ii is a measure of the efficiency of 
collision. The numerical factor in the second expressions assumes N a = 3 x 10 8 and a = lfi m. In 
this equation, n s is the density in the shell of the radiative shock ploughing through the IGM. 
Compare the above time-scale with the cooling time behind the shock (Shull & Silk 1981), 

t c ~6xl0 12 4 w s, (31) 

(1 + 8.3 x 10 3 CT~ 1/2 )n s 

where ( = Z/Z & . Therefore, dust can survive sputtering in the hot postshock region before it 
cools. 

The Lyman-a clouds are exposed to UV background radiation which can potentially destroy 
the dust particles from sublimation. However, the intensity of background radiation is several 
orders below the ambient radiation in local HII regions and it can be shown that the grain 
temperature is close to that of the microwave background radiation (Loeb & Haiman 1997; Ferrara 
et al. 1998), which is too low for sublimation to be important. 

If Lyman-a clouds are assumed to be star-forming regions then the UV flux inside these 
clouds can greatly exceed the background value (for dust heating we do not consider this case). 
From deep, serendipitous searches for redshifted Lyman-a and Ha emission lines, it is possible 
to put an upper limit on the ionizing intensity (or equivalently on the average star-formation 
rate) inside the clouds (Hogan & Weymann 1987). These searches give upper limits which suggest 
that the Hydrogen-ionizing flux inside the clouds do not exceed the background value by more 

than a factor of 10 (Thomson and Djorgovski 1995, Thomson et al. 1996). As the dust grain 

1/5 

temperature, oc JJ , this additional UV flux does not appreciably enhance it. 

In the local ISM, other destruction mechanisms are cosmic ray erosion and shattering from 
dust collisions . The properties of cosmic ray at high redshifts are largely unknown. However, even 
if we assume the cosmic ray density in the clouds to be comparable to local ISM, the timescales 
of dust erosion are > 10 10 yrs (Draine & Salpeter 1979). The timescale of dust-dust collision in 
shocks in local ISM is ~ 10 8 yrs (Jones et al. 1996). In Lyman-a clouds this time scale would 
be larger by a factor of the square of the ratio of ISM to the cloud number density of grains. As 
this factor is > 10 7 , dust-dust collision can be neglected as an important destruction mechanism. 
Also, since the magnetic field in the IGM is thought be small (B < 10~ 9 G; Kronberg 1994), 
other destruction processes like non-thermal sputtering from betatron acceleration also become 
unimportant. 
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Therefore, we assume the dust to gas ratio to be proportional to that in Galactic ISM with a 
factor C = Z/Zq. 

4.2. Lya absorbers 

It is clear from the Figures 5, 6 and 7 that dust heating can play an important role in 
determining the physical state of the gas in Lya absorbers, if they can be considered as overdense 
regions in the IGM. As mentioned in §3, the gas in Lya absorbers are most likely to behave 
somewhere between the two extreme cases considered here, expand in two dimensions while being 
confined in one direction. In this case the gas would behave close to the evolution in Model I, and 
dust can heat up the gas instead of letting it cool. 

The gas temperature in Lya absorbers, however, is not directly observable. What is instead 
observed is the line width (the b value), and lines can be both thermally and velocity broadened. 
However, Haehnelt, Rauch & Steinmetz (1996) have argued that the minimum value of b in 
its distribution corresponds to the gas temperature, at least for lines with N HI ^ 10 14 cm -2 . 
Numerical simulations seem to show that lines with lower column density are mostly velocity 
broadened (Miralda-Escude & Rees 1996). In this context, it is interesting to note the recent 
observations of Kim et al.(1997). They used Keck spectra of 5 quasars at different redshifts 
and found that the minimum value of b increased in time. We show the corresponding values 
of temperature in Figure 7 as thick lines. Although it is difficult to draw any firm conclusion, 
the evolution of temperature does suggest that dust heating can contribute towards raising the 
gas temperature in Lya absorbers as required by these observations. Note, that an alternative 
explanation for increase of the minimal Doppler width with decreasing z has been suggested 
recently by Giallongo (1997) as a result of decrease of the softness parameter in framework of the 
photoionized heating model (Ferrara & Giallongo, 1996). 

A note on the uncertainties in our calculations due to our assumptions concerning the 
background radiation is in order here. Unlike simple photoionization heating, heating due to dust 
grains is proportional to the intensity of the background radiation. This makes our calculations 
sensitive to the uncertainties in the background radiation, which include uncertainties in the QSO 
luminosity function, the contribution from young galaxies and the opacity due to Lya systems. 
As an example, using the LA model of absorption from Lya clouds of Meiksin & Madau (1993) 
(see, Sethi & Nath 1997), instead of model A3 of Giroux & Shapiro (1996), would increase the 
background flux and the dust heating rate. 

Since dust heating rate is proportional to the flux of photons, we would expect a rise in 
temperature in absorbing systems close to quasars. The uncertainties in various parameter, 
however, may make this prediction difficult to test. Another prediction of our calculation is that, 
with a large scatter in the metallicity and also in the value of Zi n , there would be a scatter in the 
temperature and in the minimum width in the distribution of b at a given redshift. 



-19- 



As we have shown in §3, if there is a low density 'intercloud' part of the IGM, then the thermal 
evolution of this component of IGM would be much different with the inclusion of dust heating 
than if only photoionization heating were used. This heating would also depend on the metallicity, 
and the metallicity of the IGM is not expected to be homogeneous. Nath and Trentham (1997) 
have estimated the probability distribution function of metallicity in the IGM, in the case of IGM 
enrichment from galactic winds, and found the standard deviation of the metallicity of order 
O.OIZq. Rauch et al.(1997) found from the comparison of data with their numerical simulation a 
scatter of an order of magnitude in the metallicity of Lya absorbers. It, therefore, is reasonable 
to assume that there are regions in the IGM with higher dust content than we have considered 
above, and so, hotter in general. (Of course, the actual temperature would depend on the thermal 
history of the region.) 

Finally, we estimate the metallicity that is needed for dust heating to be more important than 
photoionization of a primordial gas. From fig. 6, we estimate this lower limit for the accepted QSO 
type spectrum to be Z/Z & £ 10 -25 (n/10~ 5 cm -3 ) ( J9i2,-2i/0.5)~ a7 , when the grain mixture is 
dominated by silicates (as in LMC or SMC type of mixture). 



4.3. Feedback from IGM with dust heating 

The feedback of the intergalactic medium on the process of galaxy formation is an interesting 
phenomenon and has been discussed in by several authors. The reionization and reheating of the 
intergalactic medium inhibit the formation of small galaxies by the suppression of cooling and 
because of the high temperature of the IGM gas (Efstathiou 1992, Thoul & Weinberg 1996). Here 
we briefly discuss the implications of dust photo-electric heating for such feedback mechanisms. 

First, unlike photoionization heating of a primordial gas, the dust photo-electric heating is 
proportional to the UV background radiation intensity. Fluctuations in the intensity and spectrum 
of the background radiation would result in variation in the inhibiting efficiency. However, even 
for silicates and for QSO type spectrum, the gas temperature is expected to change only by 
a factor ~ 1.5, thereby increasing the Jeans mass by a factor of ~ 2. The dependence of the 
dust heating on the dust-to-gas ratio, however, has more interesting effects. When dust heating 
dominates in gravitationally confined (i.e. with constant or growing density) clouds, and the 
cooling is determined mainly by bremsstrahlung emission, A(T) oc T 1 / 2 , the gas temperature can 



be estimated from Eq. 28 as 



T OC D(«+0+l)/(a+0) j2/(a+/?) n -2/(a+/3) } ( 32 ) 

where T>, is the dust-to-gas ratio. Therefore, the enrichment of the gas (higher T>) results in 
the increase of gas temperature (and Jeans mass, Mj oc p 3 ( a +^+ 1 )/ 2 ( Q +' 3 )) ) and thus suppresses 
subsequent formation of low mass galaxies. One can speculate from these simplified arguments 
that dust heating tends to suppress the progressive enrichment of the IGM with heavy elements. 
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However, if the metallicity of enriched IGM gas exceeds a critical value ( > ( c , when radiative 
cooling from heavy elements dominates over bremsstrahlung cooling, metal-enriched regions cool 
efficiently and dust heating does not work anymore as a negative feedback. For cooling due to 
C IV, which is the most dominant heavy element for cooling in these conditions, we estimate that 
Cc ~ 0.05 (n/10 -4 ) -1 , for Jgi2,-2i = 1, a = 1, and T = 5 x 10 4 K. Note also, that dust heating can 
work as a positive feedback (stimulating formation of regions of low temperature due to thermal 
instability) if hard spectrum (a < 0.5) and large (a > 0.01 fim) silicate grains dominate, though 
this case seems to be unrealistic. 



5. Conclusion 

We summarise our findings: 

(a) It is reasonable to associate dust grains with the metallicity observed in high redshift Lya 
absorbing systems, as the probability of the survival of dust grains during the process of spreading 
the metals and dust in the IGM and in the Lya systems is high. 

(b) With the presence of the UV background radiation at high redshift, the dust particles in 
the IGM acquire charges which are much larger than is typical in the ISM, especially because of 
the low density of gas and high intensity of hard UV and X-ray photons. We illustrate this effect 
using two models of the UV background radiation, described in Table 1. Figures 1 and 2 show 
the dependences of the grain potential on gas density and grain size for graphites and silicates. 
Silicates in general have larger charges, and harder spectra help to produce larger grain potentials. 

(c) This results in a higher mean energy of the photo-ejected electrons, and therefore, yields 
an important source of heating of the IGM gas. We show the volume heating rate as functions 
of the gas density, grain size distribution and grain type. Silicates fare better as heating agents, 
and smaller grains dominates the heating rate in general. We estimate that dust heating would be 
important for our QSO type spectra for a metallicity Z/Zq ^ 10~ 2 ' 5 (n/10~ 5 cm~ 3 ) Jgi2'-2i- 

(d) Using these heating rates, and grain mixtures of the LMC type, we compute the 
temperature evolution of the gas in overdense regions of the IGM, with overdensity of ~ 10, for 
various initial redshifts, initial temperatures and metallicities. We assume that the Lya systems 
evolve in a passive manner, driven by the evolution of the background UV radiation and the 
boundary condition of the cloud expansion. We show the evolution for overdense regions which 
(i) adiabatically expand and keep the overdensity constant and (ii) which are confined either by 
shocks or gravity and keep the (proper) density constant. The difference in the temperatures in 
the cases when dust heating is included and when it is not, is of order of a factor of ~ 2 for our 
hard QSO type spectra, for a metallicity of Z/Zq = 0.01 in a Hubble time at z ~ 3 for Qf, = 0.03. 

(e) We have speculated on the connection between the dust heating in such systems and the 
observed increase in the minimum width of Lya absorption lines in time (most probably indicative 
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of increase in the gas temperature). We have also discussed the possible increase in the gas 
temperature of Lya absorption systems near QSOs due to increased UV flux, and the patchiness 
in the temperature of the IGM due to (i) density inhomogeneity and the (ii) inhomogeneity in the 
metallicity. 
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Figure captions 

Figure 1(a) - Grain potential eU/kT is plotted against ambient density n e for QSO and QSO+YG 
spectra for graphites for two values of Jgi2,-2i = 0.5, 0.05 and for a = 1CT 5 cm. 

Figure 1(b) - Same as in Fig. la but for silicates. 

Figure 2(a) - The grain potential is plotted against a for graphites for two incident photon spectra 
(QSO and QSO+YG) for J 912 _ 2 i = 0.5,0.05, and for n e = 1(T 5 , 10~ 4 cm" 3 . 

Figure 2(b) - Same as in Fig. 2a but for silicates. 

Figure 3 - Mean energy of photo-emitted electrons is shown against ambient density for a = 10 -5 
cm, for QSO and QSO+YG spectra, for graphites. 

Figure 4 - Same as in Fig. 3 but for silicates. 

Figure 5 - The mean electron energy E(a) is plotted against eU/kT for various sizes, photon 
spectra and grain types, at T = 10 4 K. The solid lines are for silicates and the dotted lines are 
for graphites. The thicker lines refer to QSO and thin lines denote QSO+YG type spectra. The 
upper set of curves are for a = 10 _5 cm and the lower set for a = 10~ 6 cm. 

Figure 6 - Integrated (over size distribution) heating rate for graphites and silicates vs ambient 
density for QSO (solid lines) and QSO+YG (dotted lines) spectra. The upper set of curves are for 
^912,-21 = 0.5 and the lower set for Jgi2,-2i = 0.05. Photoionization heating rate for a primordial 
gas is shown for QSO (short dashed lines) and QSO+YG (long dashed lines). 

Figure 7 - The heating rate is shown as a function of the lower limit of the grain size (a), for 
silicates (solid lines) and graphites (dotted lines) for two different densities n = 10~ 4 , 10~ 5 cm~ 3 , 
for QSO spectrum and at T = 10 4 K. 

Figure 8(a) - Temperature evolution for clouds with constant density (left panels) and constant 
overdensity (right panels) for graphites. The evolution is shown for Zi n = 3 & 4. The top and 
bottom panels are for QSO and QSO+YG spectra respectively. The dotted lines show the case 
with only photoionization heating, the solid lines correspond to the case where dust heating is 
included for a metallicity of Z/Z & = 0.01, for = 0.03, and Tj n = 10 4 K. Thermal history with 
the QSO spectrum is added also with Tj„ = 2 x 10 4 K. 

Figure 8(b) - The same as Fig8a but for silicates. 
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Figure 9 - Temperature evolution for clouds with constant density (left panel) and constant 
overdensity (right panel) for QSO spectrum for a mixture of grains of LMC type. The evolution 
is shown for zi n = 3 & 4. The dotted, dashed and solid lines correspond to the cases with 
only photoionization, with dust with metallicities Z/Z & = 0.003 and 0.01. The three sets of 
curves correspond to initial temperatures Tj n = 1,2,5 x 10 4 K. The thick solid lines show the 
temperatures corresponding to the minimum b values from Kim et al.(1997). 
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